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ABSTRACT 

We present the fifth release of the UMIST Database for Astrochemistry (UDfA). The new reaction network contains 6173 gas-phase 
reactions, involving 467 species, 47 of which are new to this release. We have updated rate coefficients across all reaction types. 
We have included 1171 new anion reactions and updated and reviewed all photorates. In addition to the usual reaction network, we 
also now include, for download, state-specific deuterated rate coefficients, deuterium exchange reactions and a list of surface binding 
energies for many neutral species. Where possible, we have referenced the original source of all new and existing data. We have 
tested the main reaction network using a dark cloud model and a carbon-rich circumstellar envelope model. We present and briefly 
discuss the results of these models. 
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(-H 1. Introduction 

^p^Chemical models are an important tool in helping us understand 

Q various physical and chemical processes in space. The need for 
accurate models of chemical evolution of astrophysical environ- 

c/) ments is of ever-increasing importance as a new generation of 
^ ^ ^ ground-based and space-borne facilities opens up spectroscopic 
windows at high spatial and spectral resolution. The concurrent 
development of improved receivers and laboratory spectroscopy 

J> has led to the identification of more than 150 molecular species 
CNj (ignoring isotopologues) and to the realisation that a full un- 
derstanding of the physics and chemistry in molecular sources 
requires a detailed understanding of chemical kinetics and, in 
particular, reaction rate coefficients over a wide range of tem- 
peratures: from 10 K or less in pre-stellar cores, to several hun- 
T— I dred Kelvin in hot molecular cores, to several thousand Kelvin 
in post-shock gas. 

, Although it has been widely recognised over the past decade 
J> that grain-suiface chemistry can play a significant role in molec- 

" '~j ular synthesis, the many uncertainties associated with this has 
prevented the development of accurate quantitative models for 
^ Surface chemistry. Thus, we do not attempt to include surface 

- - -chemistry, and focus on describing gas-phase chemistry as accu- 
rately as possible. 

In this paper, we present the fifth release of the UMIST 
Database for A strochemistry, Rate12 (previous releases: Rate91 
- iMillar et aL| [l991: Rate95 - Millar et al. 1997a; Rate99 - Le 
Teuff et al. 2000; Rate06 - Woodall et al. 2007). We have un- 
dertaken a major revision of the database. The gas-phase chem- 
istry is now described by 6173 reactions among 467 species, of 
which, 47 of these are new additions, composed of 13 elements. 
We have a new website (www.udfa.net ) from which the reaction 
network and associated data and codes can be downloaded. We 



* All codes, along with reaction networks and data files, are accessi- 
ble at www.udfa.net, 



also include state-specific rate coefficients for deuterium frac- 
tionation of H3 , a complete, singly-deuterated version of Rate12 
and surface binding energies of neutral species. We have made 
updates to data across all reaction types and, in particular, have 
sought out original references to the data, giving the DOI (Digi- 
tal Object Identifier) code, where available. 

We use the entire gas-phase reaction network to model the 
chemistry in two environments: a dark cloud and a carbon-rich 
circumstellar envelope (CSE) surrounding an AGB star Our 
results show that the data is sufficiently comprehensive to use 
in a range of astrophysical environments without the need to 
omit specific reactions or indeed to include additional reactions. 
Other environments for which Rate 12 can be used without modi- 
fication, along with careful treatment of molecular hydrogen and 
CO self shielding, are photodissociation regions (PDRs), hydro- 
dynamic shock regions and diffuse clouds. Protoplanetary disks 
and hot molecular cores can be modelled with the addition of 
grain-surface chemistry. 

In Sect. |2] we briefly describe the format of the data and in 
Sect. |3] we list the updates we have made. In Sect. |4] we out- 
line some of the new features of Rate12, while in Sect. |5] we 
present the results of our model calculations which we also com- 
pare with observations of TMC-1 and IRC -1-10216. 



2. Species and Related data 

We present the interstellar chemistry of 467 species, including 
268 cations, 28 anions and 171 neutral species, in the form of 
the Rate12 reaction network, which is available electronically 
as a colon-separated file. The format of each line is explained 
in Table[T] Where available, Digital Object Identifiers (DOI) are 
used for the reference codes (REF). In cases where a DOI has not 
been found or does not exist, we adopt the Rate06 referencing 
method allocating a 4 digit code to each source. The reaction 
type codes are listed in Table |2] 
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Table 1. Format of the colon-separated reaction network file 



line format: 

reaction no. :type:Rl:R2:Pl:P2:P3:P4:NE: [a:/3: y:T;:T„ : ST:ACC:REF] 



Field Heading 


Description 


Comments 


reaction no. 


Reaction number 




type 


Reaction type code 


See Table|2jfor codes 


RX 


Reactants 


Either a chemical species or one 
of: PHOTON, CRP, CRPHOT, e- 


PX 


Reaction products 


Either a chemical species or e- 


NE 


Number of fitted 
temperature ranges 


See text 


a, /3, and y 


Parameters used to calculate 
the rate coefficient 


See Sect. 12. 11 




Lower/upper temperature 


The lower/upper bound at which 
the measured or calculated data 
has been fitted. 


ST 


Source type 


One of: 

E - Estimated 

M - Measured 

C - Calculated 

L - A combination of a 

number of experimental 
values from the Literature 


ACC 


Accuracy 


Codes representing an error of: 
A - <25% 
B - <50% 
C - within a factor of 2 
D - within an order of 

magnitude 
E - highly uncertain 


REF 


Reference 


DOl or other, see text. 



Table 2. The code, reaction type and the number of each reaction type 
in Rate 12 



Code 


Reaction type 


Count 


AD 


Associative Detachment 


132 


CD 


CoUisional Dissociation 


14 


CE 


Charge Exchange 


579 


CP 


Cosmic-Ray Proton (CRP) 


11 


CR 


Cosmic-Ray Photon (CRPHOT) 


249 


DR 


Dissociative Recombination 


531 


IN 


Ion-Neutral 


2589 


MN 


Mutual Neutralisation 


981 


NN 


Neutral-Neutral 


619 


PH 


Photoprocess 


336 


RA 


Radiative Association 


92 


REA 


Radiative Electron Attachment 


24 


RR 


Radiative Recombination 


16 



For reactions where it has not been possible to define a single 
formula to fit the available data, the 'NE' field gives the number 
of different temperature ranges for which we supply a different 
formula fitting the data. Care has been taken not to have any dis- 
continuities in the rate coefficient between temperature ranges. 
Thus, if a particular rate coefficient is best fitted using two sep- 
arate Arrhenius expressions for the temperature ranges 10 - 100 
K and 100 - 1000 K, both expressions will give the same value at 
100 K. In order to evaluate the rate coefficient outside the given 
temperature range, we recommend that the user chooses the ex- 



pression that is closest to the temperature of interest. While there 
is no guarantee that this will give the correct rate coefficient, we 
have taken care to ensure that, when evaluated at low temper- 
atures (< 50 K), no rate coefficient will become unphysically 
large. We discuss this issue further in Sect. |37 



2.1. Calculation of the reaction rate coefficient using a, /3, 
andy 

For two-body reactions the rate coefficient is given by the usual 
Arrhenius-type formula 



k - a\ exp — cm^ s 

\300/ ^\ T I 



(1) 



where T (K) is the gas temperature. For direct cosmic-ray ioni- 
sation (type - CP), 

k^a s"', (2) 
whereas for cosmic -ray-induced photoreactions (type = CR), 



4—] — 

\300/ \-t 



(3) 



where a is the cosmic-ray ionisation rate, y is the efficiency of 
the cosmic-rav ionisation event as defined in equation 8 of Gre- 
del et al. (1989), and co is the dust-grain albedo in the far ul- 
traviolet (typically 0.4-0.6 at 150nm). The cosmic-ray ionisa- 
tion rates listed here are normalised to a total rate for electron 
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production from cosmic-ray ionisation (primarily from H2 and 
He in dark clouds) of ^0 = 1.36 x 10 '^i ' . Rates for both 
direct cosmic -ray ionisation and cosmic-ray-induced photoreac- 
tions can be scaled to other choices of the ionisation rate, ^, by 
multiplying the appropriate rate coefficients by ^/^o. For inter- 
stellar photoreactions (type - PH), the rate coefficient is param- 
eterised as, 

k-aexp(-yAv) s'^, (4) 

where a represents the rate coefficient in the unshielded interstel- 
lar ultraviolet radiation field. Ay is the dust extinction at visible 
wavelengths and y is the parameter used to take into account the 
increased dust extinction at ultraviolet wavelengths. 

3. Updates since Rate06 

3.1. Anion reactions 

In 2006, Ce H" was discovered in t he Taurus Molecular Cloud 1 
(TMC-n bv FMcCarthv et all (l2006tV Since then, it and other an- 
ions have_been_detectedjn et aL 
2007t ICordiner et al.ll20TTl iRemiian et al.ll2007i: iSakai et at 
2010(). Rate12 includes 22 new anions, which are involved in 
1280 reac tions. The anion reaction network comes from Walsh 
et al. (20091? with rate coefficients updated to include more recent 
laboratory measurements. 

Anions are formed primarily via radiative electron attach- 
ment, 

X -I- ^ X" -H hv, (5) 

and destruction mechanisms include mutual neutralization with 
abundant cations, photodetachment reactions and reactions with 
atoms: H, C, O and N. Of the new anion reactions, 950 are 'mu- 
tual neutralisation' reactions (A" + — » A -n B). This new re- 
action type, labelled 'MN', replaces the Rate06 ion-molecule 
and ion-ion neutralisations reaction types (IM and II). We have 
included reactions of anions with the 20 most abundant cations 
in a dark cloud chemical model, along with the top 10 most abun- 
dant cations in the CSEs surrounding both an oxygen-rich and a 
carbon-rich AGB star These reactions are all assumed to have 
the same rate coefficient, k ^ 7.51 x 10~^(r/30 0)'°^ cm^ s"', 
as included in the model of I Walsh et all (l2009l) . The products 
are obtained, where appropriate, by simple electron transfer For 
cations with no direct neutral equivalent e.g., H^, we assume the 
neutralisation of the cation results in the same products as its 
dissociative recombination e.g. CgH" H- H3 ^ CgH -H H- H 

3.2. Dissociative recombination reactions 

Since the release of Rate06, many new dissociative recombi- 
nation (DR) rate coefficients have been accurately measured, 
along with branching ratios, at the CRYRING storage ring ex- 
periment in Stockholm. In Rate12, we include new DR data for 
18 species, resulting in the addition of 58 new reactions. This 
brings the total DR product channels to 532 (compared with 486 
in Rate06), 146 of these are measured. 

3.3. Neutral-neutral reactions 

There are 74 new or updated neutral-neutral reactions in Rate 12 
which have been measured in the laboratory. These include rate 



coefficient measurements down to 24 K using the CRESU (Cini- 
tique de Reaction en Ecoulement Supersonique Uniforme or Re- 
action Kinetics in Uniform Supersonic Flow) method, and mea- 
surements of the reactions of manv hvdrocarbons (e.s. Berteloite 
et al. 2010; C anosa et ah ,2007; ,Cartv et al.i,200U Loison & 
Bergeat 2009). 

3.4. Cation-neutral reactions 

There are 174 new, and 30 updated, ion-neutral (IN) reactions 
in Rate12. The majority of the new cation-neutral reactions in- 
clude neutral species new to Rate12. In Rate06 we included 
a separate 'dipole-enhanced' reaction network, which enhanced 
rate coefficients at low temperatures in cases where the neutral 
species has a large, permanent dipole moment. Such enhance- 
ments, due to 'dipole locking' effects (Herbst & Leung 1986; 
lTYo3ri 987). result in rate coefficients which have a T '^^ depen- 
dence at low temperatures. In our reaction file we included this 
power-law behaviour in all reactions for which (i) the neutral has 
a dipole moment in excess of 0.9 Debye, (ii) the reaction does 
not already have a temperature dependence, and (iii) the reaction 
does not have a measu red rate coefficient at low temperatures. 
IWoon & Herbsj (l2009h present an alternative approach, based 
on the Su-Chesnavich expression (Su & Chesnavich 1982). In 
Rate12, we offer only the 'dipole-enhanced' reaction network 
which has been shown to give a much better fit to observations 
at low temperatures while not significantly influencing results at 
high temperature (< 1000 K). 

3.5. Radiative association 

There are just 17 updated or new radiative association (RA) rate 
coefficients in Rate12. Although these are not many in num- 
ber, several important reactions have been revisited, measured 
or calculated. These include reactions of CjHj, CHj and 
with molecular hydrodgen and the radiative association of H and 
H+. 

3.6. Cosmic-ray-induced photoreactions 

Rate12 includes at least one cosmic -ray photoreaction rate coef- 
ficient for all neutral and anionic species in the reaction network. 

3.7. Photoprocesses 

Since the last release of the database, many new photodissocia- 
tion and photoionisation rates, along with depth-dependent rates, 
have been calculated (van Hemert & van Dishoeck 2008). These 
rates along with some corrections and unpublished data, have 
been incorporated into the database (van Dishoeck, private com- 
munication 2011). We have also recalculat ed anion photo rates 
using the Draine interstellar radiation field (Draine 1978), and 
the power law described in Millar et al. (2007). Rate12 now in- 
cludes at least one photoprocess for every neutral species in the 
reaction network. 

3.8. Refits to problematic reaction rates 

In Rate06, there was a list of reaction ID's for which it was rec- 
ommended that the rate coefficients be set to zero at 10 K. This is 
because the best fit for the available data had a negative value for 
y, which, in some cases, led to u nrealistic divergent behaviour at 
low temperatures. iRoUid (1201 ih addressed this problem and re- 
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fitted many of the problematic reactions in Rate06. He fixed the 
value of the rate coefficient at 10 K to a particular value and 
used a fitting algorithm to find a new Arrhenius-type expres- 
sion for the reaction. That paper also gave new fits for reactions 
which have discontinuities at certain temperatures, or tempera- 
ture ranges in which no rate coefficient is defined in Rate06. He 
did this because, when using chemical models, it is often prefer- 
able to have a reaction rate that is continuous within the errors, 
than one that has discontinuit ies. In RATEl2 we incorporate many 
of the changes sugg ested by RoUigl (1201 ih . and we use a similar 
method to fit some of our new reaction rate coefficients. 

3.9. HNCO isomers 

Three of the four isomers of isocyanic acid (HNCO) have been 
detected in space (Briinken et al. 2010; Marcelino et al. 2008). 
A comprehe nsive gas and grain reaction set was compiled and 



modelled by lOuan et al 



gas-phase reactions from 



(120101). W e have included all of the 



Ouan et all (12010) in Rate12 



4. What's new in Rate12 
4. 1 . Deuterium cliemistry 

In order to accurately model deuterium chemistry, all D-bearing 
analogues of H-bearing species need to be included in a chemical 
model. This can increase the number of reactions by roughly an 
order of magnitude as well as increasing the number of ODEs to 
be solved. In order to reduce the number of reactions, specific 
rules can be adopted to Unut the types of reaction that can occur 
(lAlbertsson et al.ll201lHRoberts & Millajl200d ). In addition, 
branching ratios and rates can differ significantly from their H 
equivalent and are often unmeasured. For these reasons, we are 
providing a list of important deuterium exchange reactions as 
well as a singly-deuterated version of Rate12 for download from 

our website. 

Flower et al. ([20041) have shown that the H2 ortho/para ratio 
has a strong influence on the deuterium fractionation in interstel- 
lar clouds. We also include, for download, a set of state-specific 
deuterated reactions involving isotopologues of H^, H|, and Hj" , 
most ly taken fr om Flo wer et all (|2004|) . IWalmslev et al] (l2004h 
and .Hugo etalj (.2009.) . 



4.2. Reaction notes files 

We have recorded detailed sets of notes for each new and up- 
dated reaction, giving comments on possible fits and reasons for 
particular choices regarding rate coefficients. These are available 
on the website. 



4.3. Data source references 

Where possible, we have cited the original source of data for 
each reaction. Correct citing of ffie data is important, not only 
to gauge the accuracy or reliability of a particular reaction rate 
coefficient, but also to ensure the original authors are correctly 
cited in future publications. All new and updated reactions are 
referenced with a DOI (Digital Object Identifier), where avail- 
able, which allows website users to link directly to the paper's 
webpage if they require further information on the data. We have 
also provided, DOI codes, where available, for existing sources 
and located references for over 1200 reactions which were lack- 
ing source information in earlier versions of the database. In 



AppendixlAl we list separately all original sources of laboratory 
and literature data. 



4.4. INCHI codes 

We have assigned all species in the database with an IUR\C 
International Chemical Identifiers (InChI). The InChI is a lay- 
ered, variable length ASCII identifier which allows each species 
to be uniquely described. The Standard versions of InChI and 
InChlKey have been generated, where possible, for species in 
Rate12 and many of these have been cross-referenced against 
InChI data in other databases. For more information about In- 



Chls see www.inchi-trust.orgj 



4.5. VA/WDC coiiaboration 

UDfA is a part of the VAMDC project. VAMDC, the Virtual 
Atomic and Molecular Data Centre, aims to build a reliable, 
open, flexible interface to existing atomic and molecu lar (A&M) 
data hosted in various databases worldwide (Dube met et al.l 
I2OIO: Rix on et al1l201ll) . It will provide the wide community 
of both European and global users with access to a comprehen- 
sive, federated set of A&M data and application resources. For 
more information about VAMDC and its associated standards, 
see www.vamdc.eu 



4.6. Surface binding energies 

In order to model grain-surface chemistry accurately, many fac- 
tors need to be taken into account; the grain-surface morphology, 
composition, la yering and grain size all effect how the chem- 
istry progresses (ICuppen & He rbst 2005). In addition, under cer- 
tain conditions, normal rate equation methods for modelling sur- 
face chemistry are inadequate and stochastic approaches such as 
Monte-Carlo, modified rate-equation and master equation meth- 
ods are needed to mode l g rain-surface chemistry mo re accu- 
rately (iBiham et alilMTtlCaseni et al.iri998l:IChamleviri998h . 

Surface binding-energies can be used for determining ther- 
mal and cosmic -ray-induced desorption rates from interstellar 
dust-grains and also for e stimating t he diffusion barriers between 
grain surface sites (.Hasegawa et al.i,1992) . In Rate12, we have 
included a list of surface binding energies for 206 species. The 
binding energies in the list are a mixture of theoretical work, es- 
timates and measurements (Allen & Robinson 1977; Garrod & 
Herbst 2006; Hasegawa & Herbst 1993; Hasegawa e t alJl 19921 : 
iTielens & Anamandola.1987.) . .Whittet et al.. " cl998.) observed 
that water ice is the main constituent of interstellar dust- grain 
ice mantles and so, where possible, we have listed the values 
measured in water ic e. We include recen t measurements from 
iCollings et al. (2004}). IOberg et al.1 ( l2009bl) . IOberg et al.1 (l2009cl) 
and .Edridge.(,201Q) . 



4. 7. New website 



Our website Iwww.udfa.n etl has been significantly updated for 
this release of the database. In addition to being able to down- 
load the entire Rate12 reaction network, users can also now 
download the control models presented here, together with in- 
structions on how to run and modify them. A singly-deuterated 
version of the whole reaction network, a state-specific deuterium 
chemistry and a list of surface binding energies are also available 
to download. 
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Comprehensive searching of the database by chemi- 
cal species is, of course, possible, however, we have 
now integrated content with that previously available on 
|www.astrochemistry.net This means that the entry for each 
species contains considerably more information than just those 
of reactions involving it. We now also present the control model 
data for both the dark cloud and circumstellar envelope models 
and the user can also explore the effects of switching on or off 
specific reactions. In addition, it is now possible to view results 
from a sensitivity analysis on a speci es by spec i es bas is. This 
analysis follows the method of Wake lam et al.l (1201 Ol) . where 
many models are run with the reaction rate coefficients varied 
within the quoted errors using a log-normal distribution. The 
analysis provides information on the variation of the mean value 
of the fractional abundance of each species with time, and the 
first and second standard deviations away from the mean, for 
every species in the database and for both control models. 

For each reaction, as much information as possible is pre- 
sented regarding the evolution of rate coefficient data between 
previous releases of UDfA and the this one. In most cases where 
a rate coefficient has changed between Rate06 and Rate 12, there 
is now an explanatory note. In the future, changes will not be 
made without these notes, and the version of the database at the 
time of publication will remain available, even if it is not the 
current one. 

Finally, users of the website can comment directly on any 
species or reaction, for example, to discuss the values we have 
adopted or to alert us and the community to new laboratory or 
theoretical data. 



5. Results 

In this section, we use the complete Rate12 reaction network in 
two models: a dark cloud, and a circumstellar envelope (CSE) 
of a mass-losing, carbon-rich AGB star Although we shall com- 
pare our results with observations towards two archetypical and 
well-studied objects, TMC-1 (CP) and IRC-i-10216, our purpose 
here is not to reproduce the observations in detail, but to check 
(i) whether our global rate file gives a reasonable first approxi- 
mation to the chemistry in these objects and (ii) whether there 
are particular species for which our new gas-phase model dra- 
matically fails to reproduce observations. The latter item will 
highlight particular species or reactions which require further 
investigation or whether additional chemical processes are re- 
quired, e.g., grain-surface chemistry. 

5. 1 . Dark cloud model 

We model the gas-phase chemistry of a dark cloud by treating it 
as an homogeneous, isotropic cloud with constant physical pa- 
rameters: n(H2) = 10"* cm""*, T=10 K, Ay = 10 mag, a dust-grain 
albedo in the far ultraviolet of 0.6 and a cosmic-ray ionisation 
rate of 1 .3 x 10"' s"' . No grain-surface chemistry is included in 
this model except the formation of H2 via the association of two 
H atoms. This occurs at a rate of (5.2 x 10-'^)(r/300)'^-%H n(H) 
cm"^ s"' which assumes that all H atoms that stick to a grain 
surface will recombine to form H2. Our initial elemental abun- 
dances, which are listed in Table [3l are identical to those in Gar- 
rod et al. (20 08), these are the low-metal abundances of Graedel 
et al. (1 19821) . updated with more recent diffuse cloud values for 
He, C+, N and O. 



Table 3. Initial abundances relative to total H nuclei, nn 



Species / 


■ 1 

Uj/nii 


Species / 


ni/nn 


H2 


0.5 


Na 


2.0(-09) 


H 


5.0(-05) 


Mg 


7.0(-09) 


He 


0.09 


Si 


8.0(-09) 


C 


1.4(-04) 


P 


3.0(-09) 


N 


7.5(-05) 


S 


8.0(-08) 





3.2(-04) 


CI 


4.0(-09) 


F 


2.0(-08) 


Fe 


3.0(-09) 


Notes. a(b) 


= fl X 10* 







5.1 .1 . Model results 

Although we are not constructing a model of a particular dark 
cloud, it is instructive nonetheless, to compare a simple model 
with observations of TMC-1 (CP), the cyanopolyyne peak in the 
TMC-1 ridge, a dark cloud region in which over 60 molecules 
have been detected. In order to compare our model results with 
TMC-1 (CP) observed values, we used a simple method where 
we select the best-fit time as the time at which the most mod- 
elled abundances 'match' those observed. We define a match for 
a particular species when its calculated abundance agrees within 
an order of magnitude of the observed value. This definition of 
a match, while simple, has the advantage of giving equal weight 
to all species. Thus, species which are known to have an incor- 
rect or incomplete chemistry do not skew the agreement. For 
example, if we used the 'distance' method of Wak elam et al] 
(12006), the calculation of the 'best time' would be completely 
dominated by CH3OH, whose modelled abundance is more than 
3 orders of magnitude lower than its observed value at all times 
in the model. It is well known that gas-phase models do not re- 
produce the methanol abundances seen in dark clouds, so this 
method would not give a good estimate of how well the model 
fits observations. 

As our model does not include grain-surface chemistry, we 
do not expect those species thought to form in, and on, grain 
mantles, such as methanol, to have an accurately modelled abun- 
dance. Our models are also clearly unphysical for times greater 
than about 1-2 million years since freeze out onto grain surfaces 
will dominate the evolution of the gas-phase abundances. Ta- 
ble |4] shows the observed and computed best-time abundances 
for 63 species at our 'best-fit time' (defined above), ~ 2 x 10^ 
years, when 41 of 63 species have abundances falling within 
an order of magnitude of those observed. A further 14 species 
match at some time between 10"* and 10** years. In the table, the 
species listed in bold font are those for which the modelled abun- 
dance never falls within an order of magnitude of the observed 
value. Some of these species, for example, methanol and ac- 
etaldehyde, are generally thought to be produced mainly through 
grain-surface processes. However, other species, for instance, 
C2, C2H and C4H, are not thought to be influenced greatly by 
grain chemistry. The fact that these three species fall signifi- 
cantly below their observed values and the abundances calcu- 
lated with Rate06, indicates that the updated Rate12 chemistry 
results in reduced abundances of C2-bearing species and some 
of the higher-order carbon chains, whose abundances depend 
strongly on the C2 chemistry. 

One reaction which makes a significant difference to how 
well our results agree with observations, is C2H H- O — » CO + 
CH, one of the 'key reactions' in the KIDA database (Kinetic 
Database for Astrochemistry, http://kida. obs.u-bordeauxl.frj 
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Fig. 1. A plot of C2H fractional abundance, relative to H2 abundance, as a function of time using different values of the rate coefficient for the 
reaction, C2H + O — > CO + CH. The shaded region corresponds to an order of magnitude error on the observed abundance in TMC-1 (CP). 



IWakelam et al.l (l2012h ). Rate12 adopts the same high value for 
the rate coefficient, 10"'" cm^ s"', as recommended in the KIDA 
datasheets; the Rate06 value for this react ion is 1.7 x 10~" cm ^ 
s at 10 K, on the recommendation of Baulc h et a 
based on laboratory measu rements at room temperature and 
higher Subsequentlv. lBaulchl (l2005h recommended a rate coef- 
ficient about an order of magnitude larger with a constant value 
over a wide temperature range, based on a single measurement 
at 600K. As it turns out, using the 'old' lower value gives a bet- 
ter agreement with observations for a number of species such as 
C2H, C4H, CH3CHCH2 and SO2, while setting the rate coefli- 
cient to zero makes the agreement even better, this can be seen 
for C2H in Fig. [T] While one cannot accurately infer a rate 
coefficient from astrophysical data, this example highlights the 
importance of knowing rate coefficients accurately at low tem- 
peratures. 

In order to see how Rate 12 compares with the KIDA 
database, we ran a dark cloud model, using the reaction network 
'kida.uva.20ir, with the same physical conditions and initial 
abundances as above. This reaction network is the most recent 
reaction network from the KIDA website for use in dense inter- 
stellar clouds. 

As expected, the model results are largely similar to Rate12 . 
Using the same criteria as above, we find that of the 62 modelled 
species that have been observed in TMCl (CP), 38 agree at a 
'best fit' time of 2.5 x 10^ years. In Table|5] we list those species 
for which we find either Rate 12 only or kida.uva.201 1 only 
agree with observations. Of the species for which Rate12 shows 
a better fit, those that stand out include acetaldehyde, propene 
and formic acid. The differences between the kida.uva.201 1 and 
Rate12 results for acetaldehyde and formic acid can be explained 
by the newly measured dissociative recombination rate coeffi- 



cients for the reactions 

CH3CHOH+ -He ^ CH3CHO + H, and (6) 
HCOOH^ -I- e ^ HCOOH + H. (7) 

These rate coefficients are larger than their previous values and 
are each the main formation mechanism for thei r respective 
products (i Hamberg et~al]|20I0tlVigren et al.ll2010ah . Similarly, 
kida.uva.201 1 is missing reactions that synthesise propene effi- 
ciently. In Rate12, propene foiTnation is initiated by the disso- 
ciative recombination of CjH:^ as outlined in Herbst et al. (2010). 
One species that is highly underabundant in RateI2 is C2S, 
caused by the large rate coefficient of 1 x 10~'° cm~^s~' for its 
destruction by atomic oxygen. This react ion has been d eemed 
likely to proceed without a barrier bv lLoison et al.l(l2012h . and is 
another example of how lack of low temperature data can have a 
large impact on calculated abundances. 

The major difference between Rate 1 2 an d Rate06 is the in- 
clusion of anion reactions. W alsh et al.l (1200 9) investigated the 
effect of anions on the chemistry in dark clouds. They used the 
Rate06 reaction network and added a set of carbon-chain anion 
reactions. They found that the addition of these anion reactions 
resulted in enhanced abundances of several families of carbon- 
chain molecules and better agreement with the observed values 
of cyanopolyyne molecules in TMC-1 (CP). Many of the same 
reactions and species are included in Rate 12 and, as expected, 
the main differences between Rate06 and Rate12 are the same as 
those found bv Walsh et alJ (|2009). Other differences not previ- 
ously discussed are an improved agreement with abundances for 
both propyne and OCS. The former may be attributed to Rate 12 
having a more comprehensive treatment of allene (CH2CCH2) 
and propyne (CH3CCH) than Rate06, with the latter due to the 
inclusion of t he new ly reviewed and updated OCS chemistry of 
iLoisonet all (1201 2h . 
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Table 4. Observed fractional abundances, relative to H2, in TMC-1 (CP) and corresponding modelled 'best fit' values. The best fit time is 1.7 x 10^ 
years. Of the 63 species below, 41 agree, to within an order of magnitude, with observations. The 'Ref.' column gives references for the observed 
values. 
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References. (1) See'Smith et al.' ('2004') and'Walsh et al.' (12009') for source; (2)'Bell et al.' ('1999'); (3)'Kalenskii et al.' ('2004'); (4)'Suut arinen et al.l 
(2011); (5) Sakai et al. (2010); (6) Fosse et al. (2001); (7) Liaue et al. (2006); (8) Lovas et al. (2006); (9) Remiian e t al. (2006); (10) Marcelino 
et al. (2008); (1 D iSnvder et al .1 ( 12009) ; (12) lGuelin et al.l ( ll998l) ; (13) iMarcelino et alj ( 120071) ; (14) lMinh et alj 819931) 



5.2. Carbon-rich circumstellar envelope model 

In this section, we report the results of a simple chemical model 
for the CSE of a carbon-rich AGB star, using the new Rate12 
reaction network to compute the chemical evolution. For a con- 
stant mass-loss rate, the star has a uniform , spherically symmet- 
ric outflow which can be easily modelled (lMillailll998l) . 

5.2.1. CSE model details 

The model for the CSE is based on that of iMillar et all (l2000l) us- 
ing parameters appropriate to IRC -1- 10216 but ignoring the pres- 
ence of dust and gas shells in the outflow, the inclusion of which 
gives a significant improvement in agreement between models 
and o bservation in both colum n densities and spatial distribu- 
tions (ICordiner & MiIIarll2009h . The cool, expanding outer en- 
velope of the carbon-rich AGB star, IRC -1- 1 02 1 6, is one of the 
most important extra-terrestrial environments for the study of 
gas-phase chemical kinetics due to the abundance and diversity 
of the different molecu les detected in the envelope (e.g. Cer- 
nicharo et al. i2000.) . IRC -1-10216 has therefore been the subject 



of intensive ongoing astrochemical study (e.g. lAg undez et al.l 
[2012: Bieging & Rieu 1988: Cordiner & Millar 2009|.De Beck 
et al. 2012^iGuelin et al.iJ999 ). with over 80 different molecules 
detected in this source to date. We assume the c entral star loses 
mass at a uniform rate of 1.5x10"^ Moyr ' dPe Beck et al.l 
I2OI2I) . The ejected matter expands in a spherically-symmetric 
outflow at a radial velocity of 14.5 km s resulting in a gas den- 
sity distribution that falls as 1/r^, where r is the distance from 
the central star The adopted temperature profile in the envelope 
is based on an empirical fit to the gas k inetic temperature pro- 
file derived by ' Crosa s & Me ntenI (11997), the form of which is 
given in Equation 1 of Cordiner & Millar (2009). Parent species 
(with abundances taken from Agundez et al. 2012 and Cordiner 
& Millar 2009; see Table |6]l, are injected at the inner radius, 
r, = 2 X 10'^ cm, where the molecular hydrogen number den- 
sity is = 3.2 X 10^ cm"^, and the kinetic temper ature of 
the ga s is 221 K. The standard interstellar radiation field (iDraind 
Il978l) is assumed to impinge on the outside of the envelope from 
all directions, and is attenuated by a radial visual extinction of 
Ay = 6.9 mag at r,-. The radiation flux inside the envelope is 
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Table 5. A list of species whose observed abundances are matched by only one of Rate 1 2 or kida.uva.20 1 1 . The abundances shown are those seen 
the 'best-fit' time of each reaction network. These are 1.7 x 10' years for Rate 12 and 2.5 x 10' years for kida.uva.201 1. 
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calculated assuming purely absorbing grains. Extinction of the 
incident radiation is derived using the Bohlin et al. ( 1978) stan- 
dard inte rstellar gas-to-dust ratio and the extinction curve tab- 
ulated bv lSavage & MathisI (119791') . Photodissociation of CO is 
assumed to occur through absorption of 100 nm photons, using a 
single-band ap proximation to account for self-shielding ( Morris 
& Jura 1983). 

The chemical kinetic equations are solved as a function of 
radius as material traverses the CSE, until it reaches the final ra- 
dius, rf -1 X 10'^ cm, at which point the density has decreased 
to - 26 cm"^, the radial extinction is Ay - 0.02 mag and the 
majority of molecules (apart from self-shielded H2) are dissoci- 
ated. 

Table 6. Initial abundances of parent species relative to H2 



Species 


Abundance 


Species 


Abundance ' 


He 


l.O(-l) 


CH4 


3.5(-6) 


NH3 


2.0(-6) 


H2O 


1.0(-7) 


HF 


8.0(-9) 


Mg 


1.0(-5) 


C2H2 


8.0(-5) 


HCN 


2.0(-5) 


N2 


2.0(-4) 


C2H4 


2.0(-8) 


CO 


6.0(-4) 


SiH4 


2.2(-7) 


H2S 


4.0(-9) 


HCl 


1.0(-7) 


cs 


7.0(-7) 


HCP 


2.5(-8) 


SiO 


1.8(-7) 


SiC2 


2.0(-7) 


SiS 


1.3(-6) 







Notes. <" a(b) = a X 10* 



5.2.2. CSE model results 

In Table I2] we list calculated and observed column densities for 
those species detected in IRC+10216. We obtain a match be- 
tween observed and calculated column densities (to within an 
order of magnitude, as denoted by the 'Agree' column) for 31 
out of 46 of those species. In general, agreement is very good for 
the hydrocarbons, nitriles and anions, but is poor for CH3CCH 
and for the phosphorus-bearing species, which indicates that our 



understanding of the chemistry of these species in carbon-rich 
CSEs is not yet complete. 

Compared with the previous m odels for large molecules in 
IRC+10216 bv 'Mi llar et aP (l2000l) (MHBOO) and Cordiner & 
Millar (2009) (CM09), the new Rate 12 chemical network results 
in significantly better agreement between calculated and ob- 
served column densities of the larger cyanopolyynes (HC2„+iN; 
n = 2, 3, 4). The physical model and initial conditions for the 
CSE in the present model are more similar to those of CM09 
than MHBOO, so comparison with the results of CM09 is more 
relevant. Whereas, the calculated column densities for HC5N, 
HC7N and HC9N from CM09 were about an order-of-magnitude 
too large, we find that using the Rate 12 model the column den- 
sities of these species are accurately reproduced, within about a 
factor of two. The marked differences between the cyanopolyyne 
abundances calculated by the two chemical networks are shown 
in Fig.|2](top left panel), where the dotted lines show the results 
of the CM09 network and the solid lines are using Rate12 (in this 
figure, the same physical model and initial conditions are used 
for both networks). Using the CM09 network, the cyanopolyyne 
abundances peak closer to the star and reach higher peak val- 
ues. This is because CM09/MHB00 include the following class 
of neutral-neutral reactions (for n - 1-11) that are absent from 
Rate 12 (except for the case of n = 1, which is included and has 
an energy barrier of 770 K), 

C2„H + HCN HC2„+iN + H. (8) 

These reactions dominate the synthesis of cyanopolyynes in 
MHBOO and CM09 whereas in the Rate 12 model, they are pro- 
duced mainly from the reaction of CN with (poly-)acetylenes, 

HC2„H + CN ^ HC2„+iN + H. (9) 

For clarity, HC3N is not shown in Fig.|2]as it exhibits the oppo- 
site behaviour to the larger cyanopolyynes; its production rate in 
both models is dominated by Reaction |9l whereas in CM09, a 
larger photodissociation rate for the molecule, as well as the in- 
clusion of additional destruction reactions with polyynes, leads 
to a reduction in the overall HC3N abundance. 

Other species with improved matches, when compared with 
previous models, include C3N and C3O, which are greater by 
factors of 20 and 200, respectively, using the new network. This 
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Fig. 2. Top left: A plot of the fractional abundances, relative to H2, of cyanopolyynes as a function of envelope radius using the Rate12 model 
(solid lines) compared with the results from CM09 (dotted lines). Top right: A plot of fractional abundances of polyynes as a function of envelope 
radius for the Rate12 model including anion chemistry (solid lines) and excluding anion chemistry (dashed lines). Bottom left: A plot of fractional 
abundances of polyyne anions as a function of envelope radius. Bottom right: A comparison of the fractional abundances of various cations and 
electrons, including anion chemistry (solid lines) and excluding anion chemistry (dashed lines). The CgH" fractional abundance for the 'anions 
included' model is shown, for reference, with a dot-dashed line. 



brings these into better agreement with observations (see Ta- 
ble llli. In CM09 and MHBOO, C3N was produced primarily 
from HC3N photolysis, whereas in Rate12 the following reac- 
tion dominates its synthesis, 

C2 + HCN — ^ C3N + H. (10) 

The carbon-chain oxide, C3O, is assumed to originate via the 
recombination of H3C3O+ and H2C3O+ in Rate12 and CM09. 
However, the production of H3C30^ is much more rapid in 
Rate 12 due to the inclusion of the radiative association reaction, 

CzH;^ -h CO ^ H3C3O+ + /iv, (11) 

which is absent from CM09. The resulting enhanced H3C30^ 
abundance then leads to a substantially increased C3O yield. 

The sulphuretted carbon chain, C4S, is another species for 
which the calculated abundance has changed significantly from 
CM09. The abundance is an order of magnitude less when using 
the new chemistry. In the CM09 model, C4S was synthesised 
mainly via the sulphur exchange reaction, 

C4H -H S C4S -H H, (12) 

for which the analo gous reaction involv ing C2H was measured 
in the laboratory bv Smith e t al.l (l2004l) . The unmeasured C4H 
+ S reaction is not included in RateI2, and the primary C4S 
production mechanism is via the dissociative recombination of 
HC4S+. 



5.3. Impact of anions on the CSE chemistry 

Polyyne anions have been i ncluded in models of circumstellar 
envelopes since the work of iMillar et aH (l2000l) . Produced pre- 
dominantly by radiative electron attachment to neutral polyynes, 
the polyyne anions become abundant inside the main molecu- 
lar shell at a radius, r ~ 10'^ - 10'^ cm (see Fig. |2l bottom left 
panel). C(,li is the most abundant anion in the CSE as a result of 
the large CeH abundance and its large electron attachment rate 
(the bare carbon chain anion, C^, also reaches a similar abun- 
dance). 

The presence of polyyne anions boosts the abundances of 
neutral polyynes towards larger radii (~ 10'^ cm), as shown by 
the solid lines in Fig. |2] (top right panel). The main source of 
polyynes at these radii is from photodetachment of their respec- 
tive anions (as studied in the laboratory by iBest et aT.ll20lfl) . 

C„H--H/iv^C„HH-e~. (13) 

The presence of free electrons can thus have a shielding effect 
on the polyynes because in the presence of UV radiation, anions 
tend to undergo electron detachment rather than photodissocia- 
tion. 

In molecular clouds, polyyne anions can have significant ef- 
fects on the carbon-chain chemistry as a result of their high re- 
activity (ICordiner & CharnIevll20I2l;IWalsh et al1l2009h . In con- 
trast, the presence of anions in a carbon-rich circumstellar en- 
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velope has a greater impact on the chemistry as a result of their 
effect on the ionisation balance. Anions affect the ionisation bal- 
ance because they can remove a significant fraction of free elec- 
trons from the gas. For r ^ 3x10'^ - IxlO'^ cm, the presence of 
anions reduces the free electron number density by over an order 
of magnitude. A reduction in the electron number density leads 
to reduced recombination rates, and consequently more cations 
build up in the CSE (as shown by comparison of the solid and 
dashed lines for a selection of the more abundant cations in the 
lower-right panel of Fig. |2}. Within the main molecular shell, 
the majority of cations exhibit increased abundances, enhanced 
by factors ~ 10-1000. For atomic cations, however, the opposite 
effect occurs because mutual neutralisation is more rapid than ra- 
diative recombination, such that, when the abundance of anions 
becomes sufficiently large, the abundance of atomic cations fall. 
This is illustrated for the case of in Fig. |2](lower-right panel). 
The abundance of C^^H" is also shown for reference; it exceeds 
the free electron abundance by about a factor of 4 at radial dis- 
tances (3 - 4) X 10'^ cm, at which point it is the dominant charge 
carrier in the envelope. 

The increased abundances of cations have knock-on effects 
on the chemistry, for example, in the abundances of their re- 
spective recombination products. The CjHj and CjH^ cations 
recombine to C2H (plus products), so their elevated abundances 
give rise to the enhancement in C2H seen at around (1 - 3) x 
10'^ cm in Fig.|2](top-rightpanel). Other examples include C2H3 
(produced from recombination of CjH^ ), which then reacts with 
N atoms to produce a corresponding enhancement in the CH2CN 
abundance. 

Associative detachment reactions of polyyne anions with H 
atoms boost the abundances of the polyacetylenes (HC„H) due 
to the reaction, 

C„H~-i-H^HC„H-i-e-. (14) 

However, these reactions are only efficient in the outer envelope 
where the H-atom abundance becomes sufficiently large. In the 
outer envelope the density is relatively low, so the resulting con- 
tributions to the total polyacetylene column densities are small. 

6. Summary 

In this paper, we have presented the new release of the UMIST 
database for Astrochemistry, Rate12, describing, in detail, the 
updates and new additions made. We have presented results from 
a dark cloud model and a circumstellar envelope model using 
the Rate 12 network and shown that these models give reson- 
able agreement with observations, even in the absence of grain- 
surface chemistry. 

These codes, along with sample output files and instructions 
on their usage, are available for download from our website, 
lwww.udfa.net Also available for download, are the Rate12 
reaction network and files listing surface binding energies of 
species, a state-specific deuterated reaction network, a singly- 
deuterated version of Rate 12 and a list of the important deu- 
terium exchange reactions, so that users can construct their own 
deuterated networks. 

We have updated the search facility on the website, which 
now displays all available information pertaining to a particu- 
lar reaction or species. The website includes considerably more 
information than previously, including plots of molecule abun- 
dances in each model. Users can explore the effect that changing 
each rate coefficient has on the overall results of each model. 
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Table 7. Calculated and observed column densities (cm -) for a carbon-rich CSE. Of the 46 species below, 31 agree, to within an order of 
magnitude, with observations. The 'Ref.' column gives references for the observed values. 



Species 


Calculated 


Observed 


Agree? 


Ref. 


Species 


Calculated 


Observed 


Agree? 


Ref. 


C 


1.9(16)' 


1.1(16) 


Y 


1 


C2 


4.2(15) 


7.9(14) 


Y 


2 


CzH 


9.7(15) 


3-5(15) 


Y 


1 


CN 


4.1(15) 


1.1(15) 


Y 


2 


CN 


7.0(11) 


5(12) 


Y 


3 


HCO+ 


1.2(12) 


3-4(12) 


Y 


1,23 


CH2NH 


2.3(11) 


9(12) 




21 


C3 


1.8(14) 


1(15) 


Y 


1 


C3H 


1.4(14) 


3-7(13) 


Y 


1,4 


/-C3H2 


1.1(13) 


3(12) 


Y 


24 


C-C3H2 


5.8(13) 


2(13) 


Y 


1 


CH2CN 


4.1(12) 


8.4(12) 


Y 


5 


CH3CCH 


4.9(11) 


1.6(13) 




5 


CH3CN 


5.4(12) 


6-30(12) 


Y 


1,5 


C4H 


6.5(14) 


2-9(15) 


Y 


1 


C4H 


1.3(13) 


7(11) 




6 


C3N 


5.1(14) 


2-4(14) 


Y 


1 


C3N- 


1.0(12) 


2(12) 


Y 


8 


HC3N 


4.8(14) 


1-2(15) 


Y 


1 


CH2CHCN 


3.6(10) 


5(12) 




5 


C2S 


1.4(13) 


9-15(13) 


Y 


9,10 


C5 


1.5(14) 


1(14) 


Y 


1 


C5H 


4.0(13) 


2-50(13) 


Y 


1 


C3S 


1.3(13) 


6-11(13) 


Y 


9,10 


CfiH 


5.7(14) 


7(13) 


Y 


5 


CfiH 


8.3(13) 


4(12) 




5 


C5N 


3.5(13) 


3-6(12) 


Y 


1,4 


C5N- 


7.8(12) 


3(12) 


Y 


7 


HC5N 


1.4(14) 


2-3(14) 


Y 


1 


C7H 


7.3(13) 


1-2(12) 




1,4 


CsH 


1.2(14) 


5(12) 




1 


CsH^ 


2.9(12) 


2(12) 


Y 


11 


HC7N 


4.7(13) 


1(14) 


Y 


1 


HC9N 


1.3(13) 


3(13) 


Y 


1 


CP 


2.1(12) 


1(14) 




12 


FN 


5.5(09) 


1(13) 




12 


C2P 


8.0(09) 


1(12) 




13 


SiC 


9.8(12) 


6(13) 


Y 


14 


SiC2 


1.3(15) 


2(14) 


Y 


15 


SiC3 


1.6(12) 


4(12) 


Y 


16 


SiC4 


6.8(10) 


7(12) 




17 


SiN 


2.4(12) 


4(13) 




18 


SiNC 


1.2(09) 


2(12) 




19 


C3O 


6.0(11) 


1(12) 


Y 


20 


H2CS 


3.5(11) 


1(13) 




5 


H2CO 


1.5(11) 


5(12) 




22 



Notes. <" a(b) = a X 10* 

References. (1) see references in Table 5 of 'Millar et al.'('2000'); (2)'Bakker et al.'('l997); (3)'Aeundez et al.'('2010'); (4)'Cernicharo et al.' ("2000^; 
(5) Agiindez et al. (2008); (6) Cemicharo et al. (2007); (7) Cemicharo et al. (2008) ; (8) Thaddeus et al. (2008); (9) Cernicharo et al. ( 1987); (10) 
[Bell et al. (1993); (11) Remiian et al. (2007); (12) Milam et al. (2008); (13) Halfen et al. (2008); (14) Cernicharo et al. ( 1989); (15) Thaddeus 
etal. (.1,984); (16) Apponi e t al. ( 1999); (17) Ohishi et al. (1989); (18) Turner ( 1 992J; (19),Guelin e t al. (200^; (20) ,Tenenbaum et al..c2006) ; (21) 
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